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ABSTRACT 


The  optical  reflectance  of  Bi SI  has  been  measured  with  photon  energy,  -fro, 
from  2.2  to  3.6  eV.  The  band  structures  of  BiSBr,  Bi SCI .  BiSI,  BiSeBr, 

BiSeCl  and  BiSel  have  been  calculated  by  using  the  self-consistent  pseudo¬ 
potential  method.  These  results  show  that:  (i)  they  are  indirect  gap 
semiconductors;  (ii)  the  top  of  the  valence  band  is  located  at  some  k-point 
along  17  (the  z-direction)  with  finite  curvature  (small  effective  mass); 

(iii)  the  highest  valence  band  at  T  (the  center  of  the  Brillouin  zone)  is 
only  a  few  tenths  of  an  eV  lower  than  the  evergy  of  the  top  of  the  valence 
band  and  is  flat  (large  effective  mass),  especially  in  the  directions 
perpendicular  to  z.  (iv)  in  addition  to  the  above  listed  features,  the 
minimum  of  the  conduction  band  for  BiSCl  and  BiSeCl  is  located  at  some 
k-point  along  TX,  the  x-direction  and  is  a  few  tenths  of  an  eV  lower  than 
the  lowest  conduction  band  at  T. 

These  results  suggest  that:  All  the  six  compounds,  studied  in  this  project 
win  exhlbU  the  Gjjmi  effect  if  they_arejpro£erly  dO£ed.  They  are  potential 
candidates  for  microwave  oscillators. 

The  dopings  in  BiSBr,  BiSeBr,  BiSI  and  BiSel  should  be  restricted  to  p-type. 
Effective  doping  to  exhibit  the  Gunn  effect  will  be  to  substitute  either  the 
group  IV  elements  for  Bi  or  the  group  V  elements  for  Se.  For  BiSCl  and  BiSeCl, 
both  h  and  p-type  dopings  are  possible.  The  p-type  doping  should  be  restricted 
to  the  same  elements  as  suggested  for  the  other  four  crystals.  The  best 
results  for  n-type  doping  should  be  obtained  by  substituting  the  Bi  by  the 
group  VI  elements. 

A 

The  electric  field  used  for  the  Gunn  effect  should  be  applied  along  the  z- 
direction  for^the  p-type  dopings.  It  is  necessary  to  apply  the  electric 
field  in  the  x-direction  in  order  to  exhibit  the  Gunn  effect  in  n-type 
crystals. 

Comparisons  of  the  other  available  experimental  gap  energies  with  the  theoreti¬ 
cal  results  are  made.  Discrepancy  between  the  measured  and  calculated  trend 
of  gap  values  in  I-  and  Br-  compounds  is  discussed  and  a  possible  explanation 
is  provided. 

The  bonding  properties  of  these  crystals  are  presented. 
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INTRODUCTION 


The  objective  of  this  grant  (No.  77-3390)  is  to  show  whether  BiSBr, 
BiSCl,  BiSI,  BiSeBr,  BiSeCl  and  BiSel  are  candidates  for  microwave 
oscillators. 

The  criterion  for  a  crystal  to  be  used  as  a  microwave  oscillator  is 
that  it  should  exhibit  the  Gunn  effect.  Consider  a  degenerate  n-doped 
semiconductor,  then  there  are  electrons  in  the  conduction  band  even 
at  T  =  0°K.  The  current  density  under  the  influence  of  an  externally 
applied  dc  field  is  given  by 


J  =  ere  '  ( 1 ) 

where  e  is  the  electric  field,  and  a  is  the  conductivity.  Let  p  and  t 
be  the  density  and  the  life  time  of  the  electrons,  then 


2 

a  ‘ESr 
m* 


(2) 


where  m*  is  the  effective  mass  of  the  electrons.  From  the  analogy  of 


1  1  9  t 

the  free  electron  case  —  =  7  ,  the  value  of  m*  is  inversely 

m  fi2  9? 


proportional  to  the  curvature  of  the  energy  with  respect  to  the  crystal 
momentum  of  an  electron. 


As  the  electric  field  is  increased,  not  only  does  it  accelerate  the 
electrons  but  also  it  causes  scattering  of  the  electrons.  The  scattered 
electrons  can  be  in  states  such  that  their  effective  masses  increase. 
Consequently  a  decreases  and  so  does  the  current.  The  typical  J  vs  e 
curve  is  shown  in  Fig.  1.  The  region  A  indicates  that  current  decreases 
as  the  strength  of  the  applied  field  increases.  The  sample  exhibits  a 
negative  conductivity.  This  is  the  Gunn  effect. 

If  one  recognizes  that  positive  o  represents  power  consumption,  it  is 
clear  that  negative  a  means  that  the  sample  supplies  power,  and  can  be 
used  as  an  oscillator. 


From  the  above  discussions,  it  is  apparent  that  we  need  to  find  the  band 
structure,  E(R),  for  the  electrons  in  a  semiconductor,  then  decide 
whether  it  can  be  used  as  an  oscillator.  We  believe  that  this  is  the 
first  study  of  this  kind  for  the  Bi-compounds. 
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In  Section  II,  we  discuss  the  experimental  set  up  for  the  reflectance 
measurements.  The  method  used  for  calculating  the  band  structure  will 
be  described  in  Section  III.  The  results  and  the  discussions  about 
the  basic  properties  and  applications  will  be  presented  in  Section  IV. 
Finally,  in  Section  V  concluding  remarks  wil1  be  given. 

II.  EXPERIMENTAL  ASPECT 

Reflectance  data  were  obtained  for  samples  of  the  crystals  BiSI. 
Originally,  the  crystal  samples  were  in  a  matted  form  with  small 
quantities  of  reactant  materials  filling  the  interstitial  space  between 
the  actual  crystals  of  the  material  of  interest.  The  samples  were 
washed  with  a  56%  solution  of  HI,  then  methanol,  then  acetone  several 
times  each.  The  cleaned  crystal  samples  were  then  mounted  on  an 
aluminum  sample  holder  and  placed  in  the  reflectance  chamber. 

O  ° 

The  range  of  the  scans  was  from  3350  A  to  5400  A.  The  light  source 
used  was  a  tungsten  filament  lamp.  The  light  detector  was  an  uncoated 
photomultiplier  tube  (number  98438  from  EMI  GENCOM)  with  an  S-ll 
response,  operated  between  800  and  1000  volts.  The  tube  was  electro¬ 
statically  and  magnetically  shielded  and  all  the  electronics  shared 
a  common  ground  to  reduce  noise  as  much  as  possible.  A  signal  to 
noise  ratio  of  25  was  typical.  Typical  dark  count  was  0.2  microamps. 

A  more  complete  description  of  the  reflectance  chamber  and  the  optical 
system  has  been  presented  elsewhere.'  A  schematic  diagram  is  shown 
in  Fig.  2. 

The  measurements  were  performed  at  atmospheric  pressure  since  absorption 
•in  the  visible  range  is  not  a  problem,  and  no  filter  was  used  for  the 
final  runs  as  there  was  no  discernible  difference  either  way. 

The  reflectance  data  are  the  relative  reflectance  of  the  samples.  The 
ellipsoidal  mirror  assembly  was  adjusted  to  occult  the  reference  beam 
in  such  a  way  as  to  obtain  optimal  reflectances.  If  the  full  intensity 
of  the  reference  beam  were  used,  the  reflected  signal  would  be  relatively 
too  weak,  producing  a  small  and  noisy  relative  reflectance  curve.  The 
mirror  position  was  adjusted  for  each  sample  to  maximize  the  relative 
reflectance  and  to  minimize  noise, 

O 

To  perform  a  data  scan,  the  monochr Mmator  was  set  at  3350  A,  the  incre¬ 
ment  (§A)  and  the  number  of  data  po-nts  (410)  were  entered  into  the 
reflectance  program,  which  carried  out  the  scan  automatically.  The  data 
were  then  plotted  as  a  graph  of  relative  reflectance  versus  energy  or 
wavelength. 

III.  METHOD  OF  CALCULATIONS 

The  self-consistent  pseudopotential  method  is  used  to  calculate  the  band 
structures  of  the  6  Bi-compounds.  The  procedures  involved  in  the  method 
are  summarized  schematically  in  the  block  diagram  (Fig.  3). 


The  ionic  pseudopotentials  of  Bi,  Se,  I,  Br  were  determined  by  the 
local  scheme  proposed  by  Starkloff  and  Ooannopoulos.^  The  essential 
feature  of  the  scheme  is  to  multiply  the  self-consistent  atomic 
potential 3  by  a  function  which  has  the  following  form 


f  = 


1  -  e 


-Ar 


1  +  e‘X(r-rc) 


(3) 


This  function  weakens  the  potential  for  r<rc  and  approach  1  for  r>rc. 

The  A  and  rc  are  determined  by  fitting  the  valence  state  energies  to 
the  atomic  calculations. 3  Thus,  the  resulting  ionic  pseudopotential 
(effective  potential)  gives  the  same  atomic  energies  as  in  the  full 
self-consistent  calculations. 

For  S  and  Cl ,  the  extrema  of  the  outermost  s  and  p  states  are  very 
close  together  (^0. 1  A)  as  given  in  Ref.  3.  It  is  difficult  to  use 
the  f-function,  Eq.  (3),  to  fit  simultaneously  the  s  and  p  states 
energies.  By  making  use  of  the  non-uniqueness  of  the  pseudopotential 
in  the  core  region,  we  added  a  repulsive  s-potential  of  Gaussian  form 
in  the  core  region.  The  Gaussian  potential  has  the  form: 

2 

Vs  (r)  =  Ar2e'ar  (4) 

The,  charge  densities  used  for  the  initial  construction  of  the  electron- 
electron  interaction  —  Coulomb  and  exchange  potentials  --  were  obtained 
from  earlier  empirical  pseudopotential  calculations  of  the  Sb-compounds. 

We  used  the  so  called  Xa  -  exchange  of  Slater5,  and  set  a  =  1. 

The  self-consistent  cycles  were  carried  out  at  one  special  k  point  (0.25a; 
0.25b,  0.117c)  v/here  b  is  the  longest  and  c  is  the  shortest  lattice 
constant.  The  charge  densities  associated  with  this  £  point  do  not  over¬ 
lap  with  the  ones  of  4  nearest  unit  cells6. 

To  assure  the  convergence  of  the  energy  eigenvalue  with  respect  to  the 
number  of  plane  waves,  we  used  about  380  planes  waves  as  basis  functions. 
In  addition,  there  are  another  380  plane  waves  included  via  the  Lowdin 
perturbation?. 


It  took  about  8  cycles  in  each  crystal  for  the  potential  to  attain  self- 
consistency.  The  resulting  electron-electron  interaction  and  the  ionic 
pseudopotentials  are  used  to  calculate  the  band  structure,  E  vs  k. 
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IV.  RESULTS  AND  DISCUSSION 

The  measured  reflectance  of  BiSI  for  photon  energy,  -fia),  between  2.2 
to  3.5  eV  is  shown  in  Fig.  4.  The  structures  at  2.3  and  2.6  eV  in 
the  present  results  correspond  to  the  double  structure  of  2.2  and 
2.4  eV  in  SbSI.  Tentative  identification  of  the  origins  of  the  two 
structures  will  be  presented  later. 

Let's  turn  our  attention  to  the  results  of  the  sel f-consistant 
calculations.  It  is  necessary  to  show  how  good  the  ionic  pseudo¬ 
potentials  are.  We  choose  I  as  an  example  to  show  the  comparison 
of  the  ionic  pseudopotential  and  the  full  self-consistent  (core 
included)  ionic  potential.3  .The  two  potentials  are  compared  in 
Fig.  5a.  For  r  >  3.a.u.,  perfect  agreement  of  the  two  potentials 
is  obtained.  The  pseudopotential  in  the  region  r<3,  is  considerably 
weaker  than  Vionjc.  The  pseudo  wavefunctions  of  the  s  and  the  p 
states  are  compared  with  the  full  self-consistent  results  in  Fig.  5b. 
Outside  rc,  the  pseudo  wavefunctions  agree  reasonably  with  those  of 
the  full  self-consistent  calculations.  In  the  core  region,  the 
pseudo  wavefunctions  are  smoother,  i.e.  no  wiggling,  than  the  exact 
wavefunctions.  The  ionic  pseudopotentials  of  other  elements  show 
similar  features.  So,  we  simply  plotted  the  pseudopotential s  of  Bi, 

S,  Se,  Br  and  Cl,  in  Fig.  6a-6e.  The  values  of  r?  and  X  for  Bi ,  Se, 

Br,  I  are  listed  in  Table  I  along  with  the  comparisons  of  the  fitted 
valence  states  energies  and  those  given  by  full  self-consistent  cal¬ 
culations  using  the  Herman-Skil Iman  programs.  Parameters  characteriz¬ 
ing  potentials  for  S  and  Cl  are  given  in  Table  II.  These  potentials 
will  also  be  useful  for  future  calculations  of  semiconducting  solids 
with  these  elements  as  constituents.  The  Se  results  have  been  checked 
with 'the  ones  determined  by  Joannopoulos. 8 

The  band  structures  of  the  six  compounds  are  plotted  along  the  crystal 
axes  in  Fig.  7a-7f.  Since  our  present  interests  are  those  bands  near 
the  fundamental  gap,  only  6  bands  below  and  4  bands  above  the  gap  are 
given. 

There  are  several  common  features  in  these  calculations.  They  are: 

(i)  All  of  these  crystals  are  indirect  gap  semiconductors. 

(ii)  The  top  of  the  valence  band  is  located  at  some  Tf-point,  lc0,  along 
TZ  (z-direction)  and  is  a  few  tenths  of  an  eV  above  the  highest 
valence  band  energy  at  r  (T<=o). 

-A  /\  A  A, 

(iii  The  valence  band  along  the  a  (x)  and  b  (y)  directions  have  less 
dispersion  (curvature)  than  the  ones  along  the  c  (z)  axis. 


Feature  (i)  agrees  qualitatively  with  our  earlier  empirical  calcula¬ 
tions  of  Sb-compounds.4  Most  of  the  Sb-compounds ,  however,  have  the 
top  of  the  valence  bands  at  P.  The  joint  features  of  (i)  and  (ii) 
indicate  that  the  curvature  of  the  top  of  the  valence  band  is  greater 
than  that  at  r.  Consequently,  the  effective  mass  of  the  hole  at 
r  is  larger  than  that  at  the  valence  band  edge.  Coupled  with  the 
small  energy  difference  at  ¥0  and  r,  these  crystals  should  exhibit 
the  Gunn  effect  if  they  are  p-dopped  (e.g.  S  or  Se  to  be  replaced  by 
P,  or  As).  Therefore,  we  find  that  these  six  crystals  are  potential 
candidates  for  microwave  oscillators.  Finally,  the  third  feature, 
the  difference  in  dispersion  along  the  different  axes,  is  a  mani¬ 
festation  of  the  anisotropy  in  the  crystal  structure  of  these  compounds. 

The  characteristics  of  the  conduction  bands  show  more  variation  for 
the  6  crystals  than  those  exhibited  in  the  valence  bands.  The  Br  and 
I-compounds  have  the  minimum  of  the  conduction  bands  at  F,  while  the 
two  Cl -semiconductors  exhibit  the  minimum  at  some  k-point  along  TX 
(x-direction).  Therefore,  the  Br  and  I  crystals  have  the  smallest 
direct  gap  at  P  and  the  gap  for  both  Bi SCI  and  BiSeCl  appears  along  TX. 
The  values  of  the  lattice  constants,  the  fundamental  gaps  and  the 
smallest  direct  gaps  for  these  Bi -compounds  are  summarized  in  Table 
III  accompanied  by  a  few  measured  gap  values. 

The  difference  in  the  conduction  bands  of  the  Cl-crystals  from  the 
ones  of  the  other  compounds  makes  the  possibility  of  exhibiting  the 
Gunn  effect  in  n-type  BiSCl  and  BiSeCl .  In  these  cases,  the  electric 
field  should  be  applied  along  the  x-direction. 

From  Table  III,  one  notices  that  BiSCl  and  BiSeCl  have  the  same  lattice 
constants.  This  deserves  some  explanation.  We  are  unable  to  find  the 
lattice  constants  of  BiSeCl  in  the  literature.  As  a  first  trial,  the 
lattice  constants  of  BiSel  were  used.  The  self-consistent  calculations 
show  that  BiSeCl  is  a  semimetal.  Next,  we  tried  the  lattice  constants 
of  BiSeBr.  A  similar  result  was  obtained.  This  suggests  BiSeCl  may  be 
a  semimetal.  However,  we  reject  this  possibility  based  on  the  calcu¬ 
lated  total  valence  bandwidth.  Both  results  show  that  the  width  is 
about  24  eV,  which  is  roughly  a  factor  of  2  wider  than  the  other  five 
compounds.  There  does  not  exist  any  plausible  physical  mechanism  to 
cause  the  large  width  in  the  BiSeCl  except  that  the  lattice  constants 
used  for  the  calculations  are  too  large.  Finally,  we  used  the  lattice 
constants  of  BiSCl.  Results  seem  in  line  with  the  others.  However, 
as  seen  in  Table  III,  it  has  still  the  smallest  fundamental  gap.  One 
would  expect  that  the  lattice  constants  of  BiSeCl  should  be  larger  than 
the  ones  of  BiSCl  based  simply  on  the  sizes  of  Se  and  S.  The  difficulty 
in  obtaining  a  reasonable  band  structure  of  BiSeCl  suggests  that  the 
crystal  may  not  be  stable  under  normal  conditions  (atmospheric  pressure). 
Our  calculations  correlate  the  absence  of  crystal lographic  data  for  this 
particular  crystal . 


Another  point  given  in  Table  III  is  worth  explanation.  We  see  that  the 
calculated  band  gaps  of  BiSel  and  BiSeBr  are  comparable.  Because  the 
same  experimental  group  has  reported  two  different  gap  values  of  BiSel  at 
different  times,  it  is  difficult  to  argue  that  either  of  the  experi¬ 
mental  values  are  more  accurate  than  the  theoretical  prediction.  Vie, 
then,  turn  our  attention  to  examine  the  charge  density  of  these  two 
crystals.  Intuitively,  one  would  say  that  since  Er  is  more  electro¬ 
negative  than  I,  BiSeBr  should  have  a  larger  gap  than  BiSel.  However, 
v/e  found  that  this  is  too  simple  a  picture.  In  Fig.  8a-3b  we  compare 
the  total  valence  charge  densities  of  these  two  crystals  in  a  section 
consisting  of  Bi ,  Se,  and  a  halide.  The  atoms  are  labeled  in  the 
figures.  Comparing  the  charge  densities  around  the  halides,  we  see 
that  the  Br  has  more  charge  than  the  I.  This  result  is  consistent 
with  the  intuitive  feeling.  Now,  let  us  look  at  the  charge  density 
betv/een  the  Bi-Se  bonds.  The  I-crystal  has  more  bond  charge  than  the 
other.  We  can  conclude  that  BiSeBr  is  more  ionic  between  the  Bi-Br 
bond  but  less  covalent  in  the  Bi-Se  bond,  whereas  BiSel  is  just  the 
reverse.  How  does  this  bonding  property  influence  the  value  of  the 
energy  gap?  We  examine  the  charge  densities  from  the  four  highest 
valence  bands.  In  order  to  show  more  clearly  the  features,  we  choose 
a  section  (II  in  Fig.  9)  consisting  of  the  sare  Bi  and  halides  but  the 
Se  forming  the  horizontal  bond  with  the  Bi.  In  Fig.  10a  and  10b,  the 
valence  charge  densities  of  BiSeBr  and  BiSel  in  Section  II  are  plotted. 

Both  BiSeBr  and  BiSel  exhibit  more  charge  around  the  Bi  and  Se  than 
around  the  halides.  So,  the  ionic  character  of  the  I  and  Br  should 
be  manifested  by  the  charges  associated  with  lower  energy  valence 
states.  Because  of  the  stronger  ionic  nature  of  the  Br,  the  bond 
between  the Bi  and  the  Se  is  less  covalent.  As  a  result,  the  BiSeBr 
has  a  smaller  gap  than  BiSel.  Furthermore,  from  the  charge  densities  shown 
in  Fig.  10,  we  suggest  that  the  p-doping  of  these  crystals  will  not 
be  effective  for  exhibiting  the  Gunn  effect,  if  the  halides  are 
substituted  by  group  VI  elements.  It  is  necessary  to  substitute 
either  the  Bi  by  group  IV  elements  or  the  chalcogsnides  by  group  V 
elements.  For  n-doping,  only  substitution  of  Bi  by  group  VI  elements 
will  be  effective.  The  charge  of  the  lowest  conduction  band  should  be 
associated  with  the  Bi-atom,  because  as  shown  in  Figs.  11a,  b  the  corre¬ 
sponding  charge  distributions  are  around  the  metal  ion. 

Using  the  band  structure  of  BiSI,  and  the  earlier  interpretation  of  the 
Sb-compounds,  the  structures  at  2.4  and  2.6  eV  in  the  reflectance  are 
inferred  to  be  mainly  caused  by  the  transitions  from  states  near  J'o  to 
states  near  Tg  along  rZ. 

CONCLUSION 

We  have  measured  the  reflectance  of  BiSI  and  calculated  the  hand  structure 
of  tiie  six  Bi-compounds.  The  self-consistent  psejc'oootential  method  is 
used  for  the  calculations.  The  results  show: 


A.  All  the  six  Bi-compounds  are  potential  candidates  for  microwave 
oscillators,  if  (a)  all  of  them  are  p-dopped;  (b)  the  B i SCI 
and  the  BiSeCl  are  n-dopped;  (c)  BiSeCl  cannot  be  grown  under 
atmospheric  pressure. 

B.  Smaller  lattice  constants  are  necessary  for  the  semiconducting 
behavior  of  BiSeCl.  Apparently,  there  does  not  exist  any 
crystallographic  data  of  BiSeCl.  This  is  consistent  with  the 
difficulty  in  getting  reasonable  band  structure. 

C.  The  doping  should  be  specific  in  these  crystals.  For  p-doping, 
substituting  either  the  group  IV  elements  for  Bi  or  the  group  V 
elements  for  the  chal cogenides  is  the  effective  way.  The  Gunn 
effect  is  not  expected  to  be  observed  if  one  substitutes  for 

the  halides.  For  n-type  doping,  the  only  substitution  to  exhibit 
the  Gunn  effect  is  to  replace  the  Bi  by  the  group  VI  elements. 

D.  £or  p-type  doping,  the  d.c.  field  is  to  be  applied  along  the 
z-direction,  i.e.  the  direction  of  the  needle.  For  n-type  doping, 
the  d.c.  field  should  be  applied  along  the  x-direction. 

In  addition  to  the  above  findings  for  the  microwave  oscillators,  we 
also  cite  the  following  basic  properties: 

E.  These  six  compounds  are  indirect  semiconductors. 

/V 

F.  The  bands  along  the  z-direction^are  in  general  more  dispersive 
than  those  along  the  x  and  the  y-directions. 

Gi •  There  are  two  different  types  of  bonding  in  these  crystals.  Ionic 
bonding  appears  between  Bi  and  the  halides.  These  ionic  states 
form  the  lower  energy  valence  bands.  The  bonding  between  the  Bi 
and  the  chalcogenides  is  substantial ly  of  covalent  character.  The 
states  associated  with  this  bonding  contribute  to  the  bands  near 
the  fundamental  gap. 

H.  The  value  of  the  fundamental  gap  is  determined  by  the  interplay  of 
the  ionic  and  the  covalent  bondings.  More  ionic  bonding  means 
less  charge  is  available  to  form  the  covalent  bonding.  Consequently, 
the  value  of  the  gap  is  smaller. 
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Table  Captions 


Table  I.  The  values  of  A  and  rc  for  Bi ,  Se,  Br,  I  and  the  comparison  ( 

of  the  fitted  and  the  reference  atomic  energies.  j 

i 

Table  II.  The  parameters  to  characterize  the  ionic  pseudopotentials  of  Cl 

and  S.  . 

i 


Table  III.  The  lattice  constants,  the  values  of  the  fundamental  gaps,  the 
smallest  direct  gaps  and  the  available  experimental  o^p  values. 
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Figure  Captions 

Fig.  1.  Typical  J  vs.  e  curve  for  the  Gunn  effect. 

Fig.  2.  Schematic  diagram  of  system  for  reflectance  measurements. 

Fig.  3.  The  block  diagram  of  the  self-consistent  pseudopotential  method. 

Fig.  4.  Reflectance  of  BiSI. 

Fig.  5a.  Comparison  of  the  ionic  pseudopotential  and  full  self-consistent 

potential  of  I. 

5b.  Comparison  of  the  pseudo  wavefunctions  and  the  full  self- 

consistent  (exact)  wavefunctions  of  the  s  and  p  states  of  I. 

Fig.  6a.  Ionic  pseudopotentials  of  Bi. 

6b.  Ionic  pseudopotential  of  S. 

6c.  Ionic  pseudopotential  of  Se. 

6d.  Ionic  pseudopotential  of  Br. 

6e.  Ionic  pseudopotential  of  Cl. 

Fig.  7a.  The  band  structure  of  BiSBr 

7b.  The  band  structure  of  BiSeBr 

7c.  The  band  structure  of  B i SCI 

7d.  The  band  structure  of  BiSeCl 

7e.  The  band  structure  of  BiSI 

7f.  The  band  structure  of  BiSel 

Fig.  8a.  The  total  valence  charge  densities  of  BiSeBr. 

8b.  The  total  valence  charge  densities  of  BiSel. 

Fig.  9.  Atomic  configuration  in  a  unit  cell. 

Fig.  10a.  The  charge  densities  of  4  highest  valence  bands  of  BiSeBr. 

10b.  The  charge  densities  of  4  highest  valence  bands  of  BiSel. 

Fig.  11a.  The  charge  densities  of  the  lowest  conduction  band  of  BiSeBr. 

11b.  The  charge  densities  of  the  lowest  conduction  band  of  BiSel. 
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Figurc  1 
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Figure  3 


figure  6e. 
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